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Abstract 

We consider the eigenvalue equation for the Laplace-Beltrami operator acting on scalar 
functions on the non-compact Eguchi-Hanson space. The corresponding differential equa- 
tion is reducible to a confluent Heun equation with Ince symbol [0, 2, 12]. We construct 
approximations for the eigenfunctions and their asymptotic scattering phases with the help 
of the Liouville- Green approximation (WKB). Furthermore, for specific discrete eigenvalues 
obtained by a continued T-fraction we construct the solution by the Frobenius method and 
determine its scattering phase by a monodromy computation. 

1 Introduction 

The Eguchi-Hanson metric EH79] is a self-dual, asymptotically locally Euclidean (ALE) metric 
on the cotangent bundle of the 2-sphere T*S 2 . Geometrically, this corresponds to a Ricci-flat 
metric on the smooth resolution space of an A\ -singularity - such a singularity looks like the 
origin of C 2 /Z2 where the Z2-group acts by point reflection at the origin. 
Apart from Euclidean quantum gravity this metric has physical applications in string com- 
pactification. In fact, it is well known that the only nontrivial two dimensional Calabi-Yau 
manifold - this is a K3-surface - can be obtained from the Z2-orbifold limit T^/Z^ by blowing 
up its 16 A\ -singularities (where is the complex two-dimensional torus). By gluing the 
Eguchi-Hanson metric together with the flat torus metric one can explicitly construct an al- 
most Ricci-flat metric on a K3-surface Boz88 which is related to a Ricci-flat metric on K3 by 
a gauge transformation Tau82 . 

In this paper we shall examine the eigenvalue equation for the Laplace-Beltrami operator on 
the Eguchi-Hanson space. The problem is interesting from a mathematical point of view since 
the differential equation is separable and reduces to an ordinary differential equation which 
due to its singularities can be identified as confluent Heun equation with corresponding Ince 
symbol [0, 2, 12] (see |DMR78] for definitions). The Heun equation is an ordinary differential 
equation with four regular singularities on the punctured Riemann sphere. By coalescing two 
of the regular singularities to one irregular singularity one obtains the confluent Heun equation 
(analogously to the procedure by which one obtains the confluent hyper-geometric differential 
equation from the hyper-geometric one). 

The problem is also interesting from a physical point of view since the functions that are 
obtained by gluing the eigenfunctions on the Eguchi-Hanson space together with the well known 
eigenfunctions on the flat torus describe the quantum mechanical limit of string fields on K3. 
'email: malmendi@mit.edu 
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The plan of the paper is as follows: In Sect. |2]we introduce the ordinary differential equation 
that describes the radial part of the eigenvalue equation of the Laplace operator on the Eguchi- 
Hanson space. In Sect. we construct its solutions and their corresponding scattering phases 
by the Liouville-Green approximation (WKB). This extends and corrects a result in |Mig91| . In 
Mig9^ the author used the ad-hoc version of the Liouville-Green approximation. This version 
did not give approximations for the wave functions which are valid over the whole range. 
Furthermore, we give error bounds for the constructed solutions and their scattering phases. 
The explicit calculations can be found in the Appendix (cf. App. and [Bj. In Sect. E| we 
construct the exact solutions for specific discrete eigenvalues by the Frobenius method. These 
special values are given by the vanishing condition for a continued fraction. This approach is 
similar to the treatment of the generalized spheroidal wave equation in Wil28 , Lea86]. In 
addition, we can determine the exact scattering phase by a monodromy computation. Finally, 
we will show how this information can be used to compute the asymptotic scattering phase for 
this discrete set of eigenvalues. In Sect. El we present the numerical results obtained by the 
method of Sect. and Sect. and show that they match up to a high accuracy. In Sect. |H] we 
give the conclusions of this article and a brief outlook. 



2 The Eigenvalue Equation 

In this section we introduce the eigenvalue equation for the Laplace-Beltrami operator of the 
Eguchi-Hanson space. Using the SU(2) x U(l) symmetry of the Eguchi-Hanson space this 
eigenvalue equation will reduce to an ordinary differential equation of second order. A complete 
derivation of these and similar results can be found in | GM86] . |Per78j . 

The Riemannian metric on T*S 2 that is known as Eguchi-Hanson metric is Hyperkahler. More- 
over, a complex structure which is compatible with the Hyperkahler structure can be introduced 
by identifying T*S 2 with the complex manifold T*CP X . In particular, the latter can be covered 
by two coordinate charts U = U' = C 2 with coordinates ( u '>0> respectively Here, u, v! 

denote Euclidean coordinates on the base CP 1 , and £, £' parameterize the fiber of the bundle 
TCP 1 — ► CP 1 . This means that 

UDU'^C*xC, (u 1 , £') e U (1 U' : (u, £') = (u -1 , u 2 £) . (1) 

Note that T*CP X describes the minimal resolution C 2 /Z 2 of C 2 /Z,2, where Z 2 acts on (z 1 , z 2 ) G 
C 2 by (z\z 2 ) 1 ^ (-z\-z 2 ). To see this, on U := {(u,f) G U | f ± 0} use (u,£) = (^,(z 2 ) 2 ), 
and analogously for the other chart. Then (^Q) is given by z 1 <-» z 2 , and by the above we identify 

C 2 /Z 2 with T*CP X with the zero section removed. The exceptional divisor of C 2 /Z 2 therefore 
corresponds to the zero section of T*CP X . In these coordinates, i.e. outside the exceptional 
divisor, the Eguchi-Hanson metric takes the form 



EH _ Vc 4 + R 4 ( x C 4 Z lZ j 



= 



HJ R 2 \ t3 R 2 (c 4 + R i ) 

where R 2 = \z l \ 2 + |z 2 | 2 , 

and c > is the parameter of the Eguchi-Hanson metric. Introducing the Euler angles (8, (f), ip) 
on S 3 by 

z 2 J [ sinfetC^) J ' 
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we can describe coordinates (in the coordinate patch with z 2 7^ 0) by a radial part R G R>o 
and the angles 9 G (0, it), (ft G [0, 2ir), ip G [0, 4-7r). In the following, it will be more convenient 
to use r := ^/R 4 + c 4 instead of i2. 

The singularity of C 2 /Z,2 then sits at r = c. On the exceptional divisor, the Eguchi-Hanson 
metric takes the form of the Fubini-Study metric on CP 1 , and therefore assigns volume tt to 
the exceptional divisor EH79 . 
The ^-operation mentioned above takes the form 

z 1 \ ( z 1 



z 2 J " " V z 2 

In these new coordinates the metric takes the form originally found by Eguchi and Hanson. The 
Lie group 50(3) acts on itself by multiplication from the left and from the right. Let us call the 
vector fields that generate the right-multiplication (£1, £2, £3)) and the ones that generate the 
left-multiplication (£1, £21 £3)- Since 50(3) = SU(2)/Z2 = 5 3 /Z2 we have an action of these 
vector fields on the Eguchi-Hanson space. In [GM86] the authors derive that : = 

and that (£1, £21 £3j £1) define Killing vector fields that generate an su(2)l (Bu(1)r symmetry 
algebra of the Eguchi-Hanson space. 

Let A denote the Laplace-Beltrami operator on functions that is associated to the Eguchi- 
Hanson metric. Furthermore, let ^ denote a smooth function on the Eguchi-Hanson space. We 
consider the eigenvalue equation for the real positive eigenvalue 1 , i.e. (E + A)\l/ = for E > 0. 
The operator E + A can be expressed merely in terms of £ 2 = Yli=i £i > £i' i- e - 

c 4 \ d 2 /3 c 4 \ d 4£ 2 4c 4 & 



Due to their commutation relations we can diagonalize the operators £3,£3,£ 2 simultaneously 
and expand the eigenfunctions in terms of Wigner functions D\ m [9 , ift , (ft) (see |AS55j for defi- 
nitions), i.e. 



.1 



^qm 

A(j,q,(3\z) Z^n(0,^), 

j=0 q,m=-j 

where z := t-, (3 := ^-^p, and oP qm are complex coefficients, and the A(j,q,/3\z) are functions 
that depend only on the coordinate z. The Wigner functions fulfill 



£3 i m Dg m , 

£3 i q Dq m , 

where j'GN, q, m G [-j, j]n . 

Here, j, q are the quantum numbers that label the SU(2) representation. The eigenvalue equa- 
tion for \& then reduces to an ordinary differential equation for A(j, q, /3|.), 







d 2 2z d , (f3z-j(j + l))(z 
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+ -75 Z 1 + 



dz 2 Z 2 -\dz {Z 2 - l) 2 



A(j,q,P\z). (2) 



1 We are considering only the scattering case since we want to glue together the eigenfunctions on the Eguchi- 
Hanson space and the flat torus. However, eigenfunctions for E < would correspond to bounded states. 
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The differential equation © has three singularities which we have summarized in the following 
table 



z 


r 


singularity 


roots 


1 


c 


regular 


=c 2 


-1 


i c 


regular 


±4 


oo 


oo 


irregular 





(for q = the regular singularities are logarithmic). 

One should notice that our differential equation has an irregular singularity at infinity, thus is 
not of Fuchsian type. More precisely, it is a confluent Heun equation with Ince symbol [0, 2, 12]. 
To construct a continuous solution valid on the whole Eguchi-Hanson sphere we have to choose 
the regular boundary condition at z = 1 (then we can extend this solution to z = 1, i.e. the 
blown up Ai -singularity). This means that we are looking for the recessive solutions that behave 
like [z — 1)2 for z — ► 1. This means that all the eigenfunctions obtain a constant value on the 
entire CP 1 which gives the exceptional divisor of the blow-up. By approaching the singular 
point of C 2 /Z2 with different slopes, one reaches different points in the exceptional divisor. But 
for a well-defined solution of the differential equation, the limit does not depend on the chosen 
slope. 

We remark that the differential equation © does not depend on q but only on q 2 . Therefore, 
it suffices to restrict ourselves to q > 0. Then, the function \& takes the following form 

°° j j ( • ■ \ 

¥(r, 9, 1p, <f>) = Yl A & ^' P\*>) Yl [ a Qm tymtf, <t>) + vtqm ^-qnSfi, *l>, 0)) • 
3=0 9=0 m=-j 

On these functions the differential operator of Eq. (J2J is self-adjoint: Since the differential 
operator in Eq. © is already formally self-adjoint the statement follows from the application 
of the results in |CL551 Sect. 3.9] and a detailed study of the singularity at z = 1 for q = 0, 
where the differential equation is of a limit-circle type, and for q > 1, where it is of limit-point 
type - for more details see [CL55]. 

The specific problems that arise in the treatment of the differential equation @ are due to the 
structure of singularities. In particular, the major issue is to treat three singularities - one being 
an irregular singularity - at the same time. The problem is to derive the connection between 
the different bases which provide expansions of the solutions in the neighborhood of the singular 
points. Another important question is how given a system of solutions the solutions transform 
into each other when passing through a cycle around the corresponding singularity. 
For the thrice-punctured Riemann sphere with only regular singularities (this is a generalized 
hyper-geometric equation) this can be done in terms of the Meijer transcendental functions - 
as they were recently applied in |GL01j . Here, the key technique lies in a representation of the 
solutions in terms of Mellin-Barnes integrals. 

For only two singularities - one being regular, one being singular - techniques can be applied 
which are familiar from the treatment of the Bessel differential equation. This corresponds 
to the confluent hyper-geometric equation, i.e. a hyper-geometric equation where two regular 
singularities are coalescing and forming one irregular singularity. Here, the key technique is a 
generalized Borel transformation GM94b , GM94a which relates a cycle around the regular 
singularity at z = 1 to the irregular one at z = oo. 

In the case of two regular and one irregular singularities we cannot apply Mellin-Barnes inte- 
grals. Due to the irregular singularity solutions are oscillating at large real values. Thus, the 
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conditions for convergence for the Mellin-Barnes integrals when continuing into the complex 
plane to close the path of integration to a cycle are not satisfied any more. 
On the other hand, due to the third singularity at z = — 1 a cycle around infinity is also not 
homologous to a cycle around z = 1. Thus, an asymptotic expansion of the regular solution 
cannot be derived by a Borel transformation. However, as we will show in Sect. 14.21 in some 
particular cases one can still use a similar argument. 



3 WKB type solutions 



In this section we give approximations to the recessive solutions of differential equation @ that 
can be obtained by the Liouville-Green approximation (WKB). The explicit derivation of the 
results can be found in Add. IA.11 to lA.4l 

From now on we suppress the labels (3,j, q in A(fi,j, q\z). Substituting A(z) = (z 2 — 1)~ 2 w(u, z) 
in Eq. (J2J) we find 



d 



—^w(u,z) = [-u 2 f(z)+g(z)]w{u,z), 



where u 2 = (3, g{z 



dz 2 
-(z 2 - I)- 2 



(3) 



3(3+1) 



and 



(3 >»- p 
(z - a){z 2 - 1) - b 

(Z 2 - l) 2 



We remark that due to the possible values for j, q, (3 we will always have 1 < b + 1 < a and 
< a. 

Eq. (J3J) is the standard form of a differential equation of second order considered for Liouville- 
Green approximation (WKB). Following the discussion in |Ulv74j . the specific approximation 
depends on the order of the pole at z = 1 as well as the number of simple transition points 
(tps), i.e. simple zeros, in the region z > 1 which is of interest in connection with the geometry 
discussed in Sect. A simple analysis of the function / shows that we have to deal with the 
four different cases (I to IV) which are summarized in the table below. 



case 


condition 


order of pole at z = 1 


tp for z > 1 


I 


1 < b+ 1 < a 


2 


yes 


II 


b = 0,l<a 


1 


yes 


III 


b = 0,a = l 





no 


IV 


6 = 0,0<a<l 


1 


no 



For the different cases one can then apply the Liouville-Green approximation: First, one per- 
forms a Liouville-transformation of the variable z to the new variable £. The transformation is 
given by an integral equation of the form 



Vz > z Q : G(Q 



V\fW\dt 



such that C and z are analytic functions of each other, and where the function G and the point 
zq depend on the case (I to IV) we are dealing with. Simultaneously, we replace the function 

w by a function W according to w(u,z) = y\^\ W(u,£). The choices are as follows: 
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case 


G(C) 


^0 


I 


§(-0* 


> 1 such that f(zo) = 


II 




a 


III 


c 


1 


IV 


(-0* 


1 



and 



2 f a 

:=- / dtyfT®. 



The aim of this transformation is to transform the differential equation (j3J) to a differential 
equation of the form 

^Lw(u,o = (t 1 (u,o + t 2 (o)w(u,o, 

where T\ and T 2 are real-valued functions such that the approximating differential equation 
obtained by omitting T 2 has solutions which are functions of a single variable. However, this 
transformation has been done in a way that the approximate solution will still reflect the right 
behavior of the solution at the singularity at z = I and for z — > 00. From now on we suppress 
the dependence of u in w and W . 

We have found the following approximations for the solutions of differential equation (j2J that 
are regular at z = 1 



case 


approx. function 


form 


I 


Airy 


w(z) 


II 


Whittaker 


w(z) 


III 


Bessel 


w(z) 


IV 


Bessel 


w(z) 



2u(f(z) 

^-Jo(0 



'/(*) 



e l 4 Miua (2i«0 



:J0HC| 



Asymptotically, the Eguchi-Hanson metric becomes the flat metric. Therefore, the solutions of 
the differential equation © must have the following behavior 



A( z ) z ~°° \sin[2, 

Z4 



(4) 



where A, i9 is called scattering phase. Based on the results of Olver et al. (see App. IA.1I to 
I A. 41 for details), for the scattering phase we obtain the results listed below, where in the last 
column we give the equation number for the error bounds that are determined in App. 

AM 



case 



error 



I 

II 
III 
IV 



VJ3 lim 2 



-2^3^/TT^ E 
(l-a)v^ K 



2 

l+o 



+ 



+ arg T U + 



4- - 

^ 4 



1+q 
2 



2^ E 



l+a 
2 



4- - 
^ 4 



(HI 



Here, F(0, m), E(^>, m) denote elliptic integrals of the first and second kind, and K(m) = F(^, to) 
and E(m) = E(^,m) are the corresponding complete elliptic integrals (c.f. |AS55[ Sect. 17]). 
This corrects a result in |Mig91| . 
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4 Solutions related to continued fractions 



In this section we determine for which values of (3,j,q the exact solution of the differential 
equation (J2J) can be obtained by a formal power series expansion (Frobenius method) around 
the singularity at z = 1. 

For an expansion around the regular singularity at z = 1 a transformation according to £ = 
\{z — 1) is suitable. Substitution of 

A(z) = (z 2 -l)iu(() (5) 

in Eq. (J2J) yields 

= C(C + l)^(C) + (9 + l)(2C + l)n , (C) + (/3(2C + l)+/u) «(C), (6) 

where fi = q(q + l) + 1). In a neighborhood of z = 1 the two linearly independent solutions 
can be represented by the series 

oo 

U reg(0 = Ys ak( -P> j > q K k ' ( ? ) 
fc=0 

^ oo 

Wsing(C) = «r eg (C) (lnC + 5 u ) + — Y, b k{PJ,q) C k , 

^ k=0 

where B u is a real number. Since the singular solution is not unique (one can always add a 
multiple of the regular solution), the parameter B u is not uniquely determined. However, the 
parameter can be fixed by fixing the asymptotic scattering phase of the singular solution. 
By the Frobenius method (cf. ;Rab72,, Chapt. 3.6]) we obtain the coefficients (a>k)k>-i as 
solutions of the following three-term recurrence relation 

a-i = 0, 

Vk>0: a k+1 = —- — a k - — a fc _i . (8) 

(k + l)(k + q + 1) (k+l)(k + q + l) 

Notice that a rescaling of the parameter ao results in a general rescaling of all the coefficients 
( a k)k>-i since Eq. (JHJ) is linear. However, the crucial information, i.e. the ratio of a\ and ao 
is fixed by the recurrence relation l|H|). By standard methods one can show that there are only 
two types of solutions for (a k )k&n in © depending on the radius of convergence r c of the series 

lim O m = |-l r c = l, 

A solution with r c = 1 corresponds to a solution of the recurrence relation for generic values 
of ((3, j, q) whereas the solution with r c = oo is the minimal solution. 2 The question for which 
values of the parameters (/3, j, q) this minimal solution exists will be considered in Chap. 14.11 
One should mention that it is possible to derive an explicit representation of the coefficients 
(a,k)keN by the use of Babister's inhomogeneous hyper-geometric functions KxtiHa . However, 

2 A minimal solution (gn)ngN is defined by the universal property that for any other solution (/i„)„ 6 j one 
obtains lim n ^oo jf- — 0. 
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this approach did not enable us to derive the asymptotic behavior in the case of a minimal 
solution. To solve the differential equation © by a Laplace transformation as in |Ext91bj for 
similar differential equations or an Euler transformation as in |Kaz98j is not possible due to the 
structure of the coefficients in Eq. ©. 

To gain a better understanding of solutions to (JHJ), let us take a look at the solution it reg which 
is regular at z = 1 (£ = 0) in terms of the singularity at z = — 1 (£ = —1). Let i> re g(C + 1) be 
the solution regular at z = — 1 (£ = — 1). Its continuation to ( — > oo shall have an asymptotic 
scattering phase of 7, which by @ and © gives an asymptotic expansion 



WC + 1) - ^3 { 1 + °{c 



We fix a solution singular at z = — 1 = — 1) by the requirement that its asymptotic scattering 
phase differs from the regular solution by a phase of ^, giving an asymptotic expansion 



WC + 1) ~ ^3 ^ + U 

To resume, looking for the nontrivial solutions of the differential equation (J2J) we have already 
found two contributions: for a generic set of parameters (/?, j, q) the solution regular at z = 1 
(C = 0) has a regular singularity at z = —1 (C = ~~ !)• This means that it is a linear combination 
of 

^reg an d "Wsing! without loss of generality: 

Ureg(C) = cos (q) i> reg (C + 1) + sin (a) v sing (( + 1) 

C->oo 1 . /—t- . / „ / 1 



sin 




c 

Because of the singularity at z = — 1 (C = — 1) the representation of u reg by a power series 
expansion must break down for all values |£| > 1. Therefore, the asymptotic scattering phase 
of u reg can in general not be obtained from (JJJ). 

Only for those values (P n )n£N that lead to a minimal solution of the recurrence relation ([8]) 
the solution regular at z = 1 is regular at z = — 1 as well, and therefore of class C°°(R). This 
means sin a = and therefore 

Aj, ff (A») = 7 modvr. (9) 

Since we can always change the sign of the solution by an overall factor the scattering phase is 
only determined up to an integer multiple of n. In Chapt. FOl we will eventually determine 7. 
If we had made an expansion around the regular singularity at z = — 1 a transformation 
according to % = \{ z + 1) would have been suitable. Substitution of 

A(z) = (z 2 -l)U( X ) 

in Eq. (j2J then yields 

Eq. can also be understood as the isomonodromic deformation of the confluent Heun 
equation with an additional apparent regular singularity moving towards x = 1 Sla99]. By a 
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theorem Sla99 ( Sect. 2, Th.] Eq. (|10[) is then related to a Painleve equation, in our case P VI : 
The Painleve equation is the Newtonian equation of motion corresponding to the quantum 
Hamiltonian given by means of the Heun equation. 

This is clear since the 5C/(2)-invariant self-dual metrics are specified exactly by a solution of 
the Painleve equation P VI Oku03 . 

4.1 Minimal solutions and continued fractions 

According to Pincherle's Theorem (see |JT80l Sect. 5.3]) a three-term recurrence relation 

V„>1 y n+ X = SnUn+Jnyn-l 

has a minimal solution if and only if the following continued fraction converges: 

OO I 

7fc 7i 



72 



5 2+ 73 



6 3 + . . . 

In particular, if (a n ) nS N is the minimal solution it follows that ao ^ and 

OO I 

a ^ \S k 

For any set of parameters (j, q), where as before q > 0, we now define the function A4(j, q \x) 
by the continued fraction of Thron type (or T-fraction) 

OO I 

M(j,q\x) := So+J^TT' 
k=i ' k 



where Vfe> 5 k 
Vfe>i 7fc 



k(k + 2q+ 1) +/i- 
(k + l)(k + q+l) 
2x 



(k + l){k + q+l) 

Using the Umordnungssatz |Per771 Kap. 6.42, Satz 2] and JT80, Sect. 7.3, Th. 7.23], we can 
conclude that the function M(J, q\x) is a meromorphic function on C. Moreover, by the above 
its real zeros determine the values of x = — (3 for which the recurrence relation (JHJ) has a minimal 
solution, i.e. the radius of convergence of the power series expansion (JJJ) becomes infinite. 
The corresponding solution u reg is then of class C°°(R). Since the relation A((3,j, q\z) = 
A(—(3,j, q\ — z) holds for the differential equation ((21), it follows that any such smooth solution 
for (3,j,q in z is simultaneously a smooth solution for —(3,j,q in the variable —z. Therefore, 
x = (3 must be another zero of A4(j,q\x). This is 

M(j,q\P)=0 M(j,q\-f3)=0. 

Using the Umordnungssatz |Per771 Kap. 6.42, Satz 2], one can see that one can cancel the factor 
(k + l)(k + (7 + 1) in 7^, 5k- Namely, it is equivalent to calculate the zeros of the meromorphic 
function A4(J, q \ x) instead of q \ x) where M(j, q\x) is defined by the continued fraction 



00 

\ - c k \ 



M(j,q\x) := do + ^r^, (11) 

fc=i lClk 

where Vfc>o d k ■= k(k + 2q + 1) + /j, — x , 
Va:>i c fc := 2k(k + q)x. 
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For this continued fraction one can even prove separate convergence (see |Thr91| for definitions) . 
To see this, recall that [i = (q—j)(j + q+l), and rewrite the continued fraction in the following 
way: 

Vj_ 3 >i M(j,q\x) = d + TT i 2 oU ~ g) . 

k=1 \ak m(j,q\x) 
x 

m(j,j\x) 



Vj> M{j,j\x) = x\-l + 

Here, m(j, q\x) := Yl'kLj-q+i can ' 3e wr hten as a T- fraction using again the Umordnungssatz 
|Per77l Kap. 6.42, Satz 2], i.e. 

°° Fk x\ 

mu,q\x) = Yl irrfc' (12) 



where \/ k >j- q +\ G k 



, , - + G k x 

k=j-q+l 

1 



(k- j + q)(k+j + q + l) 

2k(k + g) 



Ffc : (A ._ i + g _ 1)(A ._ j + g)(A;+i + g)(A . +i + ?+1) . 
#/-<H-l := J - 9 + 1 • 

Since we have both YlT=j-q+i l-^fcl < 00 anc ^ EfcLi-o+i l^fcl < 00 we can a PP^ |Thr9H Sect. 
3, Th. 3.1] to the T-fraction (|12j) : let A k (x) and B k (x) be the numerators and denominators, 
respectively, of the n th approximant of the T-fraction (|T2*|) . Then the sequences (A k (x)) k >j~ q +i 
and (-Bfc(x))fc>j_q + i converge, uniformly on compact subsets of C, to entire functions A(x) and 
B(x) of order at most one. Further -B(O) = 1, -A(O) = 0, ^4'(0) = Fj^ q+ \ so that neither function 
is identically zero, and -m(j, q\x) is well defined at x = 0. 

4.2 The determination of the scattering phase and the monodromy 

In this section we calculate the scattering phase of the regular solution u reg in the case that it 
is also regular at £ = — 1. 

Let us first look at the asymptotic expansion of the solutions of Eq. 0. For an asymptotic 
expansion a Fabry transformation |Ulv74j is suitable, i.e. a change of the variable according to 
x 2 = £. With n(C) = U{x) Eq. © becomes 

(x 2 + l)U"(x) + ("^^ + 2(g + 1) ^ 2 + 1) ) U'{x) + (4/3(2x 2 + 1) + 4 M ) U(x) = . 

By standard methods (c.f. CL55j) one finds that this equation has two linearly independent 
solutions H\,H2, such that for some 5 > the following asymptotic expansions hold: 

V, -2, + 5 < argx < ir - 5 : Hl (x) ~ * E , 

k= °oo (13) 
V x + 5 < argx < 2vr - 5 : ff a (s) ~ E > 
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where the coefficients (pk)keN fulfill the following four-term recurrence relation 



P-2 = p-i = , 

Vfc > : = — — p fc + 32/3 — — p k _ x 

K + 1 K+l 

k 2 + (2q - l)k - q + \ 
- 32/3 k+ \ Pk-2- (14) 

It follows from this recurrence relation that for po E K the coefficients (pk)keN are real. In 
particular, we have H\(x) = H2(x). This property remains true not only asymptotically, but 
also for the actual solutions: Since all the coefficients in the differential equation are real, the 
complex conjugate of any solution is again a solution. 

The Stokes phenomenon will lead to a non-trivial monodromy of the solutions if we pass through 
a cycle around infinity Q i— ► C,e 2m . This is because for x \— > xe m we will leave the sector in which 
the asymptotic expansion ((11)) for H2 holds |DU95j . However, since the solutions H\,H2 are 
a complete system of solutions the new solutions can be expressed as linear combinations of 
them. These are the connection formulae. In the case of a second order differential equation 
they had been explicitly calculated, namely in [DO 94]. We obtain for \x\ 3> 

limHiixe™ 6 ) = -2irP H^x) + % H 2 (x) , (15) 

6 — *1 

\imH 2 (xe* i6 ) = iHUx). 
e-*i 

The parameter P can be determined by a generalized Borel transformation of the asymptotic 
solution as pointed out in |(irM94aj . |GM94bj . For ordinary differential equations of second 
order this has been done explicitly by Daalhuis and Olver |DU94j . Applying these results we 
can determine the parameter P from the coefficients in the asymptotic expansion 1)14)1 . i.e. 



fe^oo [K — 1)1 Pq 

In particular, the parameter P is real. We subsume the solutions H\,H2 in the vector H := 
\HuH2). 

We need to calculate the linear combination of the regular solution u rcg (() in terms of the 

solutions H{x) derived earlier. This amounts to determining a complex parameter A E U(l) 
such that 

^reg(C) = (A X)-H(X). 



The occurrence of A and A is due to the fact that w reg (C) is real, i.e. u reg (C) = u reg (C)> and 
H^)=H 2 {x). 

Now, the representation of u reg (C) by the power series (JJJ) is valid for |£| > 1 iff the coefficients 
(afc)fc>i constitute a minimal solution of the recurrence relation (jSJ). Therefore, w reg (C) has a 
trivial monodromy around both of the regular singularities iff the 1 constitute a minimal 

solution of Eq. (jSJ). Hence, iff for any set of parameters (f3,j,q) the equation M(j,q\f3) = 
holds, then the following equation must hold: 



» lcg (C) = lim n rcg (Ce 2 ^ e ) = ( A A ) • UmH(xe 7 ' 



(A A)-( ~ 2 f I ).H(x). 
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Thus, the equation M(j,q\f3) = holds iff 



If we set A = i e 17 with 7 E R we obtain (since P is real) 

7 = — mod 7T , P = 




(17) 



Thus, A encodes the crucial information for the asymptotic expansion of the regular solution, 
i.e. 



This is the desired formula for the scattering phase in the case that the set of parameters (J3, j, q) 
induces a minimal solution of the recurrence relation JHJ). 

4.3 Computation 

The numerical determination of the scattering phase consists of two steps: First, we have to 
determine the successive zeros of the function M(J, q\x). However, one should mention that for 
the necessary evaluation of the continued fraction one has to use a backward algorithm since 
any forward algorithm must be numerically instable as shown in |Gau67j . 
We have determined the positive zeros (P n )nen of M(j,q\x) by a simple bisection algorithm 
where the evaluation of the continued fraction was accomplished by the Gautschi algorithm. It 
is essential to use the function M.(j,q\x) instead of the earlier defined M(j,q\x): Because of 
the structure of the coefficients in the continued fraction the typical values of A4(J, q \ x) become 
very small, and the determination of its zeros unstable. With ((3 n ,j, q, 7) we then have the data 
needed for a numerical interpolation of the scattering phase over the whole range of (3. 

5 Numerical Results 

5.1 Numerical Results for the WKB approximation 

The figures below show the scattering phases for different quantum numbers j and q with 

varying parameter (5 obtained by the WKB approximation from Sect. |31 

In Fig. ^ t° Fig. we have applied the WKB approximation of case I. 

Fig. n shows (from the top to the bottom) the graphs for Aj j<? for j = 1, . . . , 10 and q = 1. 

Fig. [21 shows (from the top to the bottom) the graphs for Aj i9 for j = 3, ... ,10 and q = 3. 

Fig. 01 shows (from the top to the bottom) the graphs for Aj <q for j = A, ... ,10 and q = 4. 

Fig. ^ shows (from the top to the bottom) the graphs for Aj <q for j = 5 and q = 1,3, 5. 

In Fig. |S]and Fig. Elwe have compared the WKB approximation of case II to IV (left) with the 

WKB approximation of case I (right) for the values j = 9, q = 0. The vertical line in Fig. [5] is 

indicating the transition from case II to IV. 

Fig. 13 shows (from the top to the bottom) the graphs for Aj t o for j = 1, 3, 5. 
Fig. [BJshows the graph for Ai2,o- 
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5.2 Comparison between WKB and Probenius method 

In this section we compare the numerical results for the asymptotic scattering phase obtained 
by the WKB approximation (already shown in Sect. I5.1|) with the numerical data obtained by 
the Frobenius/continued- fraction method as described in Sect, "4*"""1 

The figures below show the scattering phases for different quantum numbers j and q with varying 
parameter j3. Each figure will show the graph (/3, A^ KB (/3)) of the corresponding WKB approxi- 
mation, and crosses will mark the points ([3 n , Aj r q ((3 n )) n >i obtained by the Frobenius/continued- 
fraction method where (3 n is the n th zero of A4(j, q\(3) and 

a£,(ao = -\- h,M * - max (°>i - 1)^ (is) 

with V n >i : lj. q (n) £ No- Notice that the scattering phase is only determined up to a multiple 
of 7r since we can always change the sign of a solution by an overall factor. Therefore, Eq. (|18l) 
is equivalent to Eq. (|17[). We will also present an additional and more significant diagram of 
^Aj j(? (/?) mod 7r^ and varying parameter [3. 

Figs. 131 fm "T~""l "T~""l show Aj t o(/3) with varying j3. For j < 6 it is l{n) := n and all /3 n > j(j + 1), 
i.e. 1 > a > and we have to compare with case IV of the WKB approximation, j = 7 is the 
lowest value where for q = we have (3\ < j(j + 1), i.e. 1 < a, and we have to compare with 
case II of the WKB approximation: we have to relabel the (3 n . This can be done by setting 
Vj>7 : l{n) := |n — || — ^. 

However, since the phase can be obtained only up to multiples of tt the function l(n) is irrelevant 
for the comparison of Aj^ KB (/3) and Aj r (/3 n ). For this purpose one can look at Figs. HUllT^lfTll 
1161 which are independent of l{n). These figures show an exact match of the data sets obtained 
by the WKB approximation and the Frobenius/continued- fraction method. 

Figs. 1171 to l2H show A J) i(/3) with varying (3. For j < 8 it is l{n) := n. j = 9 is the lowest 
value where for q = 1 we have (5\ < /3 max , where /? max denotes the local maximum of the WKB 
approximation: We have to relabel the [3 n . This can be done by setting \/j> 9 : l{n) := \n— 1| — 5. 
Figs. 1171 to OH show an exact match of the data sets obtained by the WKB approximation and 
the Frobenius/continued- fraction method. 

Figs. 051 to 13*21 show Aj ) 2(/3) with varying j3. For j < 11 it is l(n) := n. Figs. H5\ to "3*21 show an 
exact match of the data sets obtained by the WKB approximation and the Frobenius/continued- 
fraction method. 

The accuracy of the WKB approximation can be explained as follows: from the error bound in 
Eq. (|29j) we can deduce the following numerical result for the error 5j tq (P) of Aj )(J (/3) 

|<5,-, 9 =o(/3)| < I min(l,l.levf -l) . 

This means that for (5 > 4 we already have 

\S M =o(l3)\ < |(eW - l) , 

and for (3 > 100 we find \5 j>q = (P)\ < ig- 

The same computations can be made for q > 1 with similar results since we also have a similar 
error bound for the general 5j, q {(3) in Eq. (J!"""]) . 
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6 Conclusions and outlook 

The eigenvalue equation for the Laplace-Beltrami operator acting on scalar functions on the 
non-compact Eguchi-Hanson space reduces to a confluent Heun equation ((J2J) or after a suitable 
substitution ©) with Ince symbol [0, 2, 12]. 

With the help of the Liouville-Green approximation (WKB) we have constructed approxima- 
tions for the eigenfunctions by special functions in Sect.|3J Depending on the quantum numbers 
that label the S*C/(2)-representation the approximating functions are either Airy, Whittaker, or 
Bessel functions. Furthermore, we have derived the scattering phases and error bounds in these 
cases. 

Moreover, for specific discrete values of the eigenvalue in Sect. H]we have constructed the exact 
solutions by the Frobenius methods. These eigenvalues are given by the zeros of a meromorphic 
function defined by the infinite continuous fraction Hll|). Together with a monodromy compu- 
tation this has provided us with the data needed for a numerical interpolation of the scattering 
phases. 

Finally, in Sect. [5] we have shown that these two sets of data (obtained by the WKB approx- 
imation and the Frobenius/continued- fraction method) agree to a high accuracy. This shows 
that one can find a discrete set of exact values for the spectral density of the Laplace-Beltrami 
operator by the method described in Sect. 14.31 Conversely, it shows that the expressions for 
the eigenfunctions and scattering phases which were derived in Sect. 131 and which have the 
advantage of being given in terms of explicit functions are very accurate approximations and 
can be used for all numerical purposes. 

It is now interesting to ask whether the meromorphic function defined by the continued fraction 
(fTTj) can be expressed explicitly as a ratio of special functions. If so we could obtain the discrete 
eigenvalues - for which we have already calculated the exact scattering phase and which we 
have also computed numerically by the continued fraction - as zeros of this ratio of special 
functions. A first step towards an explicit representation of the continued fraction might arise 
from the method of Pincherle. If we apply |Per77| Kap. 21.84, Satz 8] to the continued fraction 
(|11|) we obtain the following result: a representation of the meromorphic function defined by 
the T-fraction (|11|) is given by 



equation shows that such a solution always exists. A further analysis of the behavior at the 

irregular singularity z = 00 then shows that all the appearing integrals also exist. 

A closer investigation of the meromorphic function q\x) is subject of our ongoing research. 
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A The different cases in WKB approximation 

In this appendix we show the explicit construction of the approximate solutions of the dif- 
ferential equation @ by the Liouville- Green approximation (WKB) in the cases I to IV (cf. 
Sect. |SJ). The construction of an error bound in cases I to III can be found in the references. 
The construction of the bound in case IV is given in App. Eland might explain the construction 
and philosophy behind these error bounds. 

To give bounds for the different cases I to IV we need the notion of auxiliary weight, modulus, 
and phase function E, M, and 9: If X,Y are solutions of the respective differential equation of 
the first and second kind, then X = M E^ 1 s'm9 and Y = EM cos 9. 

In particular, in case I we will need E, M, 9 for the Airy function (see |Qlv741 Sect. 11.2] for 
more details): 

Mix) 

Mix) = — ^ sm9(x) , Bi(x) = E(x)M(x) cos9(x) . 
E(x) 

In the cases III and IV, we will use the functions Eq, Mq, 9q for the Bessel function (see |Ulv74l 
Sect. 12.1, 12.3] for more details): 

Jo (a) = cos 9 (x) , Y (x) = E (x) M (x) sm9 (x) . 

For case II, the definition of modulus and weight function is the most complicated. Therefore, 
we refer to |Dun941 Chap. 2, Th. 1] for the quite extensive definitions in this case. 



A.l Case I 



In case I it is easy to prove that the function / has a transition point (i.e. simple zero) at z > 1 
which we denote by zq. The idea is now to perform the transformation of the variable z and 
the function w(z) to (, W(() according to 



and w{z) = 
£ — > oo to z 



Vz > zq : 
\/z < zq : 



(-c) : - 



2a§ 



L z VW)dt, 



^ W{Q. Note that ( — > —cc corresponds to z — > oo, ( — > 0— to z — > zq+, and 



1. Eq. Q becomes 



HO 



16C 2 



[4f(z)f"(z)-5f'(zf 



c 



Cg(z) 

16/(z)3 f(z) 



(19) 
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Approximate solutions of ©, i.e. solutions of (fT§)) with ip = 0, that are regular at z = 1, are 
given by |()lv74l Sect. 11.3, Th. 3.1] 



w(z) 



c 



with \e(u, C)| < 



1 M(u§C) 



Ai(n3C) + £ ( u ,C) 



(20) 



J ^v flOO (|CpB (0) _j 



The constant A and the function are defined as follows: 



A := sup<J ir\x\-M 2 (x) } , 



MO 

and V is the variational operator, i.e. 

V Ci oo flCl^o(C) 



dv 



\V\ 2 



dv 



-\i/>(u,v)\ . 



C F 2 



Following the discussion in |Ulv741 Sect. 13.7.2] we can determine the behavior for z — ► 1 and 

z — ► oo: 



A(z) 
A{z) 



Z->1 



2— >oo 



1 

3" 

7T Z4 



(z-1) 



sm u 



where we have used lim^oo e(n, Q = 0, and the phase 5 is determined by lim^-oo e(u, 0. We 
know that the argument of the sin-function equals 2y / /3z + Aj j9 for large z, i.e. 



1 



3 

7T Z* 



sin ( 2v73z + A jtQ 



Therefore, for the scattering phase we obtain 

A jiq = u Urn y/f(fy dt - + ^ + 5 , 



and with the help of |Ulv74l Sect. 11.2] it follows that 



2|*| 



7T 



< min < 1, 



A 



(21) 



To resume, the solution (j2Qj) can be considered as relation between the recessive solution and 
its asymptotic expansion. In this sense, the solution is a connection formula and it is known as 
Gans-Jeffreys formula. 
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A.2 Case II 

In case II, f(z) = Jzjj^, we can apply a result of Olver and Nestor that has been generalized in 
|Dun94j . The Liouville transformation takes the form 



t — a 



(It 



with a 



2 f a 



7T 



Therefore, we perform the transformation of the variable z and the function w(z) to £, W(£) 
according to 

i i a ( 2C-a + 2C^(C-a)^ 
C2(C " a)2 -2 ln a = 



and w(z) 
kind: 



' ' lz W{C). The integral can be expressed through an elliptic integral of the second 



2d {Z - a){Z + 1) -2VTT-aE 
z — 1 



arcsm 



z — a 
z-1 



l + a 



Note that £ — > a corresponds to z — > a, £ — * to z — > 1, where branches of £ must be chosen 
such that it is an analytic function at both values. Moreover, £ — ► oo corresponds to z — > oo. 
Eq. © becomes 



' 2»-c i , V>(0 



C 4C 

4C 2 + a 2 



c 



w(0, 



+ 



16C (a - C) 2 

(C - «) g{z) 



+ [4/(z)/"(*) - 5/ 



'^21 



16/(z)* 



/(*) 



Approximate solutions of (J3J), i.e. solutions of (|19|) with tjj = 0, that are regular at z = 1, are 
given in |Dun94| Chap. 2, case 1, m = 0] in terms of the Whittaker function M, i.e. 



w(z) = 
with \s(u, Q\ < 



Vc/w 



Ml 1 2 n (2<) 

2 ' U 

£^ n (2«C) 

9 > U 



e 4 



2tt 



M 



2 ' u 



1 + e nua 
i 

: ^V 1: , (0) (F)_ 1 



(2i<) + e(n,C) 



(22) 



For the error bound we have adopted the notation of |Dun94L Chap. 2, Th. 1]. Using formula 
(2.19) in |Dun94j we determine the behavior of the solution for z — > 1 



A(z) 



2tt 



1 + e" 



2cm 2 



a — 1 V z — 1 



c 
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where we have used lim^o e(u, Q = 0. A close examination of the Liouville transformation 
then shows that (/(z - 1) -»■ (o - l)/(2a) for ( ->■ and z ->■ 1. Using jDun941 (2.16)-(2.18)] 
we obtain for z — > oo 



A(*) 



2 

2 4 



— sin(uC-— ln(2«C) + - + 7 + <y 



where 



. 1 iua\ 



and the phase <5 is determined by lim^oo e(u, Q. We know that the argument in the sin- function 
equals 2y/~/3z + A^q for large z, i.e. 



A(z) 



— sin ( 2 V/3z + Aj )0 | a >i 



A careful examination of the Liouville transformation reveals the following relation between 
large £ and z 



C-|-|ln(C) + |ln(^) + 0(r 1 ) = 2v / i-2VTT^ E 



l + a 



+ 0(z-2) . 



Therefore, for the scattering phase we obtain 



o 



Aj, | a >i = «lim C--ln(2<)-2^ +-+7 + * 



7T 



-2u v 7 ! + a E 



, 1 iua 

arg T — | 

2 2 



l + a 



au au ,au. tt 

H ln( — ) H h7 

2 2 v 2 ; 4 



2|o| 



7T 



< min < 1, 



v „(o)0F) 



(23) 



The phase A^o | a >i differs from the one obtained in Mig91 



A.3 Case III 



In case III, the function / has neither a transition point nor a pole for z > 1. This case can be 
understood as the limit a — ► 1 of case II, and applying further results of Olver's |01v77j . in this 
limit the Whittaker function will become the Bessel function Jo- 

The idea is to perform the transformation of the variable z and the function w(z) to £, W(C) 
according to 

C = / \/W)dt = 2^z + l - 2V2 



and w(z) 
becomes 



a| W(C). Note that ( — » 00 corresponds to — » 00 and £ — > to z — ► 1. Eq. © 



<K0 



4C 2 C 



W(C) 



(24) 



1 

4C 



+ [4/(z)/"(z)-5/'(z) 2 ] 



c 



Cg(f) 

16/(*)3 ' /(*) ■ 



+ 
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Approximate solutions of (J3J), i.e. solutions of with ip = 0, that are regular at z = 1, are 
given by |Qlv771 Chap. 5, Th. 2, case 2, z/ = 0, /x = 0], i.e. 



with |e(u, C)| < 



1 



Z ,i MoK) 



(2 J (u()+e{u,() 
e lo,oV ox (H) _ 1 



(25) 



h,o E (u() 

The constants £o,0) ^0,1) an d the function H are defined as follows: 



^o,o : 
^0,1 : 

H(u,C) : 
and V is the variational operator, i.e. 



sup{tt^o(x)M 2 (x)} , 

x 

sup{7rfi (aO |Jb(a?)| E (x)M (x)} , 



«o(C) 
1 + x 



1 /" c ciu 



2 y Oo(™) 

We determine the behavior of the solutions for z — > 1 and z — > oo: 



A(z) 



2— >00 



V2 



3 

7TUZ4 



sm u 



^W)dt + - + 5 



where we have used lim^ + e { u i = 0, and the phase 8 is determined by lim^oo e(u, (). We 
know that the argument of the sin-function equals 2\fjJz + A^q \ a= \ for large z, i.e. 



A(z) 



z— *oo 



V2 



3 



sin [isffiz + A Ji0 | a =i 



Therefore, for the scattering phase we obtain 



Aj, |o=i = u lim ( / v / 7M*-2\/i) +7+ S 



and with the help of |Qlv77| Chapt. 5] and analogous to App.[B]the bound for e in (|25|) gives 
a bound for the phase 5. It follows that 



< mm < 1 , — — 

[ /o,0 



•Vo.ooCfl) _ 1 



(26) 



The phase A^o differs from the one obtained in Mig91 by j. 



19 



A.4 Case IV 



In case IV the function / has no transition point for z > 1 and takes the form f(z) = Z -T\ - The 
idea is now to perform the transformation of the variable z and the function w{z) to £, V^(C) 
according to 



(-0 



dt 



and w(z) = JZd*w(£). The integral can be expressed through elliptic integrals of the first 
and second kind: 



VfUjdt 



1 z 2 — 1 

2\l + (l-a)v / 2F 

z — a 



arcsm 



z-1 



z — a 



l + a 



2V2 E 



arcsm 



z - 1 



z — a 



l + a 



Note that ( — > —00 corresponds to z — > 00 and £ — ► 0— to z — > 1+. Eq. Q b 



ecomes 



d? 
dQ 



2^(0 



u 



1 m 



w(0, 



4C 4C 2 C 
1 o(z) 4/(z)/"(z)-5f(z) 2 



(27) 



16C 4/(z) 64/(z)3 

Approximate solutions of ©, i.e. solutions of (|2*7|) with ifi = 0, that are regular at z = 1, are 
then given by j()lv74t Sect. 12.4, Th. 4.1, v = 0, n = 0], i.e. 

1 



ICI^JoHCI 5 ) + e(«,C) 



(28) 



with |e(u,C)| < t — \(\- 



e ^v c , (|ci^o(C)) _ 1 



A ,o Eq(u() 

The bound for the error e is derived in App.[Bl The constants Ao,0)Ao,l, as well as the function 
£>o are defined as follows: 



A: 



0.0 



sup {ttxMq (x)} 



x>0 



A0.1 := sup {ttxMq(x) , cos 0q(x) } , 

x>0 



Bb(C) == 



1 



dv 



r V^) , 



and V is the variational operator, i.e. 

V c ,of|Cl^o(C) 







dv 



'C M 2 

We determine the behavior of the solutions for z — ► 1 and z — > 00: 



A(*) 
A(*) 



^ -^=J (n^2(l-a)^^T 

2— »00 1 



1 



3 

TTUZ1 



sm u 



v / 70^ + ~ + <5 
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where we have used that limf_»o- £ M C) = and the phase 5 is determined by lim^—oo e(u, C). 
We know that the argument of the sin-function equals 2y / /3z + A^q for large z, i.e. 



A(z) 



1 



3 

7TUZ 4 



sin (2yf/3z + A ii0 



Therefore, we obtain for the scattering phase 

A i|Q | a<1 = u\m^(^j^ yfj(t)dt-2^j +^+ 5 



(l-a)y/20 K 



1 + a 



2V2/3 E 



In addition, we will derive in App. \B[ 



2\S\ 



7T 



< min < 1 



M),l 
V),0 



3 ^V_ 00 ,o(|Cl^o(0) 



1 + a 



+ I + 6 



(29) 



The phase Aj >0 | a<1 differs from the one obtained in [Mig91| by j and the factor in front of the 
elliptic integral of the first kind. 



B Error Bounds 



In this chapter we apply the construction of error bounds presented in |Ulv74l Sect. 11.3] to 
case IV (cf. j()lv741 Sect. 12.4, Ex. 4.4]). In case IV the Bessel functions |C| 5 Jo MCI and 
|C| 3 YoMCI 5 ) °f fi rs t and secon d type, are first approximations for W. They are exact solutions 
if ifi = 0. Substituting W(() = |C| 1 JoMCl 5 ) + e M C)> the differential equation becomes 



d 2 

^MC) 



!L _ J_ 

4C ~ 4(2 



e(n,C) 



^(0 
C 



e MC) + ICI 5 JoMCI 5 ) 



If we rewrite it as an integral equation we obtain a Volterra integral equation 



e(u,C) 
where K((, v) 



C 



7T f 2 



W 2 



e(f) + |u| 2 J («|v|2) 



JoMCI 5 )YoM^I 5 ) - YoMCI 3 ) JoMH 5 ) 

Using the fact that Eq is a non-increasing function we obtain 

VC<«<0: |K(C,u)| < M ° (n|Cl 1 2 ' > nlv^Moiulvl^Eoiulv^) 



E (u\(\2 
--■■Po(0 



:Q(v) 



We introduce the parameters Ko, k by 
k : = 



sup {P ( W )Q( W )}<^ 
«6[C,0] U 



sup 

»e[C,o] 



Q(u) \v\z Jq(u\v\2) I 



< 



il 
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With the help of |01v74l Sect. 6.10, Th. 10.2] which establishes a bound for the solution of a 
Volterra integral equation we end up with a bound for e, 



K 

< — 



k(^C)l 

Po(C) ^ 



exp{«oV f ,o(|C|'Bb(C))}- 

^ U{^v Ci0 (|ClH(C))}-i 



A, 



o.o 



(30) 



We are interested in an error bound for the phase in the asymptotic expansion of the solution. 
Assume that the solution W takes the following form for large arguments £: 



W(u,C) c ~°° 



2ich 



7TU 



1 + p) sin(n|C|5 + 1+8) + (1) 



As in |Qlv74| Sect. 6.7] we rewrite the difference of the exact and the approximate solution as 
a trigonometric function 



e(w,C) 



^(n,C)-|CI 5 JoHCI 3 ) 



(,—>oo 



2ich 



nil 



2ICM 



7TU 



{[l + p] sin («|C|3 + j + "Sin (™lCI^+£)} 



a sin [u\C\2 + ~ + , 



where the new parameters a (always considered to be positive) and r\ are related to p and 5 by 

(l + p) e iS = l + ae ir > . 

By elementary geometry it follows that 2l£L < a . Then, the bound for e in gives a bound 

for a. Now, we choose a sequence (£ n ) with linv^oo ( n = — oo for which u|C| 5 + f + ?7 is an 
odd integer multiple of 1|. This shows that 



Finally, since 



iru 



lim 



\e(u,C)\ ■ 



we obtain a bound for 5: 

TT 



M Or) 
E (x) 



2 \ - 

7TX 



< min < 1 



M),l 



A 



0.0 



e " 



•V-cc,o KI'Bo(C) 
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Figure 1: A jtQ for j = 1, . . . , 10, q = 1 
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Figure 3: A j>q for j = 4, . . . , 10, q = 4 
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Figure 5: A 9)0 for case II, III, IV 
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Figure 7: A^o for j = 1, 3, 5 
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Figure 2: A,,, for j = 3, . . . , 10, q = 3 
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Figure 4: A,-,, for j = 10, g = 1, . . . , 10 
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Figure 6: Ag ; o for case I 
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Figure 8: Ai 2 ,o 
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Figure 17: A™ and A sv 




Figure 18: A™ an d Ajjjj 
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Figure 20: A™ and A & 
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Figure 22: A^ KB and A^ 
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Figure 24: A^ KB and Agj 
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Figure 25: A™ and ^ 
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Figure 27: A™ and A^ 2 
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Figure 29: A^ 2 KB and A^ 2 




Figure 26: A™ and Agj 
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Figure 28: A™ and A^ 2 
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Figure 30: A^ 2 KB and A^ 2 
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